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SUMMARY 

In  January  1964  an  investigation  was  begun  to  investigate  the 
validity  of  classical  rigid  aircraft  analytical  techniques  on  predicting 
the  stick  fixed  and  dynamic  longitudinal  stability  of  the  Goodyear  YAO-3G 
Inflatoplane.     In  addition,   a  study  was  made  on  the  effects  of  fuselage 
flexibility  on  the  dynamic  response  in  the  phugoid  mode  and  the  movement 
of  neutral  and  maneuver  points. 

It  was  found  that  rigid  techniques  were  valid  for  dynamic  response 
predictions  and  somewhat  less  valid  for  static  predictions.     However,    if 
the  airframe  pressure  was  maintained  at  the  designed  level  static  pre- 
dictions  were  in  no  way  adverse.     If  the  airframe  is  available  for  static 
bending  test,    correction  factors  may  be  obtained  that  when  applied  to 
rigid  values  of  Cm      result  in  more  accurate  predictions  of  both  dynamic 

h 

and  static  stability. 

Fuselage  rigidity  was  varied  by  changing  internal  airframe  pres- 
sure from  5  psig  to  7  psig.     The  effect  of  fuselage  flexibility  was  to  re- 
duce damping  and  increase  the  frequency  of  the  phugoid  mode  to  a  slight 
extent.      The  neutral  and  maneuver  points  both  moved  forward;  however, 
it  was  indicated  that  N       moved  further  forward  than  N    . 


INT  ROD  U  C  T  IP  N 

The  Goodyear  YAO-3G  is  an  aircraft  that  was  developed  in  the 
195  0's  and  is  no  ionger  in  production.  The  aircraft  is  unique  in  that  it 
is  of  inflatable  rubber -fabric  construction. 

In  October  of  1963  the  Department  of  Aerospace  and  Mechanical 
Sciences  of  Princeton  University  obtained  one  of  these  aircraft  on  loan 
from  the  U.    S.   Army.     It  was  determined  that  the  unique  features  of  the 
aircraft  could  be  put  to  good  use  m  an  investigation  of  first  order  aero- 
elastic  effects  due  to  structural  bending.      The  investigation  would  be 
simplified  by  the  low  aerodynamic  pressure  required  to  manifest  such 
bending.     In  addition,   the  rigidity  of  the  aircraft  could  be  varied,   with- 
in limits,   by  controlling  the  internal  pressure  of  the  airframe. 

With  these  possibilities  in  mind  this  investigation  was  started  in 
January  of  1964.  The  first  flight  test  was  conducted  on  January  3,  1964 
and  the  tests  continued  through  early  May  of  1964. 

The  purpose  of  the  investigation  was  twofold.     First  to  investigate 
the  validity  of  classical,    rigid  aircraft,   analytical  techniques  in  pre- 
dicting static  and  dynamic  longitudinal  stability,    and  secondly  to  investi- 
gate the  effects  of  fuselage  and  wing  bending  on  the  stick  fixed  maneuver 
and  neutral  points  and  the  dynamic  modes  of  longitudinal  motion  of  this 
aircraft. 

The  investigation  consisted  of  analytical  developments  consider- 
ing the  aircraft  first  rigid  and  then  semi-rigid  and  comparing  the 
analytical  analysis  with  flight  test  results. 


EQUIPMENT 

The  Aircraft 

The  YAO-3G  is  a  single -placed  high-wing  monoplane  of  inflatable 
fabric  construction  (Figures   1  through  6  and  Table  I).      The  landing  gear 
consists  of  a  single  fixed  wheel  with  wing  tip  and  tail  skids.      The  aircraft 
is  powered  by  a  single  engine  driving  a  fixed  pitch  propeller.      The  engine 
is  a  Nelson  Model  H63-A  which  develops  a  maximum  of  42  brake  horse- 
power at  4000  RPM  at  sea  level.      The  engine  operates  on  the  two-cycle 
principle,    has  a  single  magneto  ignition,   and  air  pump  for  controlling 
airframe  pressure. 

Other  than  being  inflatable,    the  airframe  is  of  conventional  form 
consisting  of  a  straight  constant  chord  wing  (NACA  0015),   a  fuselage  with 
partially  enclosed  cockpit,    and  empennage  with  vertical  and  horizontal 
stabilizer  and  control  surfaces.     A  conventional  stick  and  rudder  pedal 
control  system  is  used  for  flight  and  manual  control  is  provided  for  the 
engine  mounted  air  pump  and  airframe  pressure  relief  valve. 

The  rigidity  of  the  aircraft  is  maintained  by  hoop  stresses  result- 
ing from  internal  pressures  varying  from  4.  8  to  7.  0  psig.     The  shape  of 
the  non -cylindrical  parts  is  maintained  by  drop  stich  threads  running 
from  one  surface  to  the  other  as  shown  in  Figure  6.     The  internal  pres- 
sure changes  with  altitude  and  is  maintained  at  the  desired  level  by 
manipulation  of  the  pump  control  and  airframe  pressure  relief  valve. 

The  instrumentation  is  quite  basic  consisting  of  an  air  speed 
indicator,    altimeter,   accelerometer,    magnetic  compass,   air  pressure 
gage,    engine  tachometer,    and  cylinder  head  temperature  indicator 
(Figure  7).     The  electrical  system  consists  of  a  dry  cell  battery   powered 
warning  light  for  the  air  pressure  gage. 


The  basic  weight  of  the  aircraft  was  315  pounds  including  test 
instrumentation  and  parachute.     Gross  weight  varied  from  500  to 
600  pounds    depending  on  pilot  and  fuel. 

Flight  Test  Instrumentation 


The  instrumentation  requirements  were  to  provide  the  necessary- 
data  to  perform  static  and  dynamic  longitudinal  stability  analyses.  This 
required  a  method  of  measuring  and  recording  elevator  deflection,,  nor- 
mal acceleration,  airspeed,  angle  of  attack,  and  pitch  rate.  In  addition, 
a  movie  camera  would  be  desirable  to  study  fuselage  bending  during  the 
phugoid  motion,  and  voice  communications  requirements  dictated  self  - 
contained  radio  equipment. 

The  instrumentation  problems  were  somewhat  different  than  would 
be  expected  of  a  more  conventional  aircraft.     The  space  available  for  in- 
strumentation was  a  small  area  behind  the  pilot  on  the  cockpit  floor.     This 
space  had  to  contain  the  major  instrumentation  components  and  power  sup- 
ply.     Weight  was  also  a  problem  in  that  a  2  00  pound  pilot  plus  parachute 
and  nominal  fuel  put  the  aircraft  at  near  maximum  gross  weight.     In 
addition,   attachment  points  were  limited  and  some  method  had  to  be  de  - 
vised  to  attach  and  remove  instrumentation  components  without  damage 
to  the  fabric  airframe. 

A  SFIM  A20  Flight  Recorder  was  chosen  to  record  data.      This 
instrument  is  6. 5x4.  7x3. 9  inches  and  weighs  4.  6  pounds  including  air- 
speed and  acceleration  transducers.     A  matched  rate  gyro  3x3x3  inches 
weighing  3  pounds  was  used  to  record  pitch  rate.     Transducers  in  the  re- 
corder convert  inputs  to  mirror  deflections.     The  mirrors  reflect  light, 
from  a  single  source,    to  light  sensitive  paper.     Five  input  channels  may 
be  recorded  on  the  60  mm  paper.      The  instrument  and  associated  equip- 
ment draws  2  amps  at  28  volts  D.  C. 


The  A2  0  can  be  set  to  run  at  tape  speeds  of  1  mm  or  5  mm  per 
second  and  holds  approximately  25  feet  of  recording  paper,     To  compen- 
sate for  variations  in  tape  speed  with  supply  voltage,   an  internal  clock 
mechanism  is  provided  which  projects  time  marks  on  the  trace.      In  addi- 
tion,   a  trace  marking  device  is  incorporated. 

The  recorder  transducers  consisted  of  an  airspeed  indicator  that 
required  static  and  dynamic  pressure  inputs,   a  self-contained  accelero- 
meter,    two  ratiometers,   and  a  galvanometer  for  potentiometer  inputs. 
The  ratiometers  compensate  for  variations  in  supply  voltage  and  were 
used  for  pitch  rate  and  elevator  deflection  inputs,     The  galvanometer  is 
not  insensitive  to  supply  voltage  and  was  used  for  the  angle  of  attack  in- 
put in  that  a  qualitative  value  of  angle  of  attack  would  be  sufficient  for 
the  dynamic  test,     A  more  complete  description  of  the  SFIM  A20  may  be 
found  in  Reference  1. 

The  power  supply  consisted  of  24  nickel  cadmium  batteries  with  a 
nominal  output  of  29  volts  and  capacity  of  4  amp  hoars.     The  power  sup- 
ply was  made  up  in  two  packs  of  12  batteries  each  in  order  to  better 
utilize  space  requirements. 

The  mounting  problem  was  solved  by  the  use  of  Permacel  P69 
fabric  tape.     This  tape  was  found  invaluable  for  manufacturing  attach- 
ment points.     It  was  found  to  adhere  well  through  the  various  tempera- 
ture and  dynamic  pressure  conditions  encounted  during  flight  test. 

Figure  8  shows  the  installation  of  the  A20  recorder,    rate  gyro, 
and  power  supply. 

In  order  to  accurately  measure  elevator  deflection,    it  was  neces- 
sary to  put  a  potentiometer  on  the  elevator  itself.     Due  to  fuselage  bend- 
ing,   the  elevator  could  deflect  with  the  stick  fixed  and  cable  measurement 
was  not  considered  to  be  of  sufficient  accuracy.      The  potentiometer 


mounting  is  shown  in  Figure  9.      The  potentiometer  ieads  were  taped  to 
the  fuselage  right  side.     Due  to  the  low  dynamic  pressure  of  the  aircraft, 
external  mountings  were  not  a  problem. 

The  mounting  of  the  angle  of  attack  vane  presented  a  problem  as 
to  what  part  of  the  aircraft  would  be  rigid  enough  to  use  as  a  reference, 
A  typical  wing  mounted  boom  was  out  of  the  question  due  to  wing  twist 
and  flexibility.     As  a  compromise  solution  the  boom  was  mounted  by  a 
two-foot  bracket  attached  to  the  landing  gear  frame.     The  gear  frame  is 
the  most  rigid  part  of  the  aircraft;  however,    in  order  to  preserve  the 
rigidity  of  the  mounting  bracket,    the  boom  had  to  be  shorter  than  neces- 
sary to  totally  escape  from  wing  interference.     However,    it  was  con- 
sidered that  the  amount  of  wing  interference  encountered  would  not  ad- 
versely affect  angle  of  attack  dynamic  response.     Although  primarily 
designed  as  a  mounting  for  the  angle  of  attack  vane,    the  mounting  bracket 
seemed  to  be  the  ideal  location  for  the  movie  camera  {a  surplus  gun 
camera),   also.     Later  on  it  was  also  used  to  mount  a  venturi  tube.    This 
venturi  was  necessary  to  amplify  static  pressure  for  the  recorder  air- 
speed transducer.      The  transducer  was  found  to  be  insensitive  at  the 
low  airspeeds  encountered  when  attached  to  a  simple  pilot  static  tube. 
Figure  10  shows  the  angle  of  attack,,    camera,    and  venturi  mounting. 

A  control  panel  for  the  recorder  was  mounted  on  the  instrument 
panel  in  place  of  the  magnetic  compass  (Figure  7).      This  panel  contained 
a  master  switch  to  control  power  to  the  recorder  motor  and  lamp,   a 
speed  switch  that  controlled  the  recorder  motor,   an  event  marker  switch, 
and  a  five  position  selector  switch  and  meter.     The  selector  switch 
allowed  the  meter  to  read  either  elevator  position,    angle  of  attack,    pitch 
rate,   battery  voltage,    or  off.     This  selector  arrangement  proved  in- 
valuable in  that  recorder  inputs  could  be  checked  in  flight.     In  addition, 
a  switch  panel  was  mounted  on  the  right  cockpit  wall.     This  panel 


contained  angle  of  attack,    pitch  rate,   and  camera  switches.     The  position 
of  the  angle  of  attack  and  pitch  rate  switches  determined  whether  these 
values  were  being  recorded.     As  these  values  were  not  necessary  for  the 
static  test,   their  presence  made  trace  reading  more  difficult.     The  switch 
panel  allows  static  and  dynamic  tests  to  be  conducted  on  the  same  flight 
without  excessive  traces  on  the  tape.     When  the  camera  switch  was  on, 
the  recorder  motor  switch  also  controlled  the  camera.     An  event  marker 
switch  was  also  mounted  on  the  top  of  the  pilot's  control  stick.     Figure  ^P-    // 
is  a  wiring  diagram  of  the  complete  instrumentation  package  and  Figure  12 
is  a  schematic  of  the  recorder  system. 

The  voice  communication  system  consisted  of  separate  self-con- 
tained VHF  transmitter  and  receiver.     The  transmitter  was  a  Navguide 
Corporation  Model  AT -103  which  was  5x7x2  inches  and  weighed 
5  pounds  including  batteries.      The  Navguide  is  a  crystal  tuned  5  tube 
VHF  transmitter  which  operates  on  dry  cells.      The  filaments  do  not 
light  until  the  microphone  button  is  depressed.     It  was  found  that  the 
filament  battery  would  not  sustain  semi-continuous  operation  necessary 
during  flight  test.      To  aleviate  this  the  recorder  power  supply  was  put 
in  parallel  with  the  filament  battery  after  being  dropped  down  to  14  volts 
with  a  4000  ohm,    5  watt  resistor.     The  transmitter  was  mounted  behind 
the  pilot  just  under  the  leading  edge  of  the  wing.     The  transmitter  antenna 
was  taped  to  the  left  outer  cockpit  wall  (Figure  3).     The  receiver  was  a 
Montery  Electronics  "Partner"  continuous  tuning  VHF  receiver  2x4x5 
inches  weighing  If;  pounds  including  mercury  batteries.     The  receiver 
was  mounted  on  the  left  inner  cockpit  wall  (Figure  7).      The  receiver 
antenna  was  taped  to  the  left  outer  cockpit  wall  (Figure  3). 

As  the  aircraft  has  no  metal  skin  there  was  no  ground  plane  effect 
for  the  antenna  mountings  and  the  receiver  would  not  operate  under  these 


conditions.     A  ground  plane  was  constructed  from  1/16  inch  welding  rod 
and  taped  to  the  underside  of  the  cockpit.     This  arrangement  proved  quite 
satisfactory  for  two  way  communication  up  to  5  miles.     As  this  range 
was  satisfactory  for  flight  test  no  ground  plane  was  used  for  the  trans- 
mitter antenna. 

The  total  weight  of  the  instrumentation  including  all  mountings 
amounted  to  33.  5  pounds.     The  SF1M  recorder  operated  very  well 
throughout  the  flight  test.     Its  light  weight,    compact  size,   and  simple 
power  supply  requirement  solved  the  basic  instrumentation  problem. 


ANALYTICAL  PROCEDURES 

Longitudinal  Stick-fixed  Stability 

The  longitudinal  stability  is  determined  by  considering  the  motion 
of  the  aircraft  in  horizontal  and  vertical  translation  and  in  pitch.     Symmetry 
of  the  aircraft  precludes  any  coupling  between  longitudinal  and  lateral 
motions.     The  wind-axis  is  assumed  to  be  parallel  to  the  principal  longitu- 
dinal axis  in  the  vertical  plane  of  symmetry,    thus  eliminating  the  products 
of  inertia.      The  resulting  perturbation  equations  of  motion  are,    from 
Reference  2. 


— ^  AV  +  Atf   +  ~  La    -  A0  =  0 


MgA6  =  M6   A6e 

In  this  evaluation  the  aircraft  will  be  considered  as  a  rigid  body 
except  where  full-scale  wind  tunnel  test  data  is  available.      Using  the  air- 
craft specifications  of  Table  I  and  references  as  noted,   the  stability  deriva- 
tives may  be  evaluated  as  follows.     The  flight  conditions  used  for  evalua- 
tion will  be: 

W  =  515.2  lbs. 
c.  g.  =  29%  mac 
h  =    sea  level 
V  =  50  knots  =  84.  45  fps 
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The  lift  coefficient  is  determined  from: 


.      _     2L  (2)(515.2) 

JL      PV2S      '(.  00238K84.45)3{124.3) 


The  drag  coefficient  is  determined  from: 


CD  "  CD     +  eTTAR 


With  Cn     =  .  1  from  Reference  3,   and  a  reasonable  value  of  the  aircraft 
^o 

efficiency  factor,    e,    taken  as  .  75  <,   the  drag  coefficient  becomes: 
„      _     ,    x  (.488)2 


D       '  (.75)(3.14)(4.44) 


The  aircraft  non-dimensional  mass  based  on  wing  chord  is: 


=  10.28 


pSc-        (.00238)(124.3)(5.31) 


The  effect  on  the  forces  and  moments  due  to  a  change  in  forward 
speed  with  angle  of  attack  and  throttle  fixed  are  stated  below.     Drag 
changes  with  speed  by: 


11. 


pSV  _  (.00238)(124.  3)(84,45)(,  123) 

DV  "     m    °D  ~  ~  16~~        "       ~ ~ 


The  thrust  versus  velocity  curve  shown  in  Reference  3,    gives  a  change  of 
thrust  with  velocity  of  -«  711  lbs/ ft/  sec. 


V       m  dV  16 


-.  04445 


The  change  of  lift  with  a  change  of  velocity  is  expressible  by: 


V2       (84.45)3 


■—  =  .  00903 


With  the  engine  lever  arm  3.  14'  above  the  center  of  gravity  and  assuming 
the  thrust  line  parallel  to  the  fuselage  reference  line,  the  change  in  pitch- 
ing moment  with  speed  is: 


M  1    SM  _  1     BT  (-.  711)(-3.14) 

V      ly     SV  "Iy   bV  277.2 
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The  changes  that  occur  in  the  forces  and  moments  with  a  change  in 
angle  of  attack  are  normally  an  increase  in  lift,  am  increase  in  drag  and  a 
negative  pitching  moment.  The  lift  curve  slope  of  the  aircraft  at  5  0  knots 
as  determined  from  full-scale  wind  tunnel  tests,   Reference  4,    is: 

CL      -  4.18     /rad 
a 

The  drag  change  with  angle  of  attack  is  represented  by: 

D  2S        r         =     (2)(32.2)(4.  18) 

a     eTTAR  ^L  (.  75)(3.  14)(4.  44) 

D      =  26.  0 
a 

A  change  in  angle  of  attack  effects  the  lift  by: 
a  _        ^a   _    (32.2)(4.  18) 


L 

^=3.265 


In  determining  the  angle  of  attack  effect  on  the  pitching  moment  the  con- 
tributions of  the  fuselage  and  tail  must  be  considered  first.      The  fuselage 
contribution  may  be  estimated  from  Reference  5  by: 

KfW/L£ 
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Where  L.    is  the  overall  fuselage  length9   W    is  the  maximum  width  of  the 
fuselage,   and  K    is  an  emperical  factor  determined  from  experimental 
evidence. 


r  _  (.012)(2.25)3(19.6) 

Umof  "        (124.3)(5.31)  "  ' 


An  angle  of  attack  change  at  the  tail  produces  a  moment  change  as  follows: 


The  tail  efficiency,   r\    -  -*—    is  felt  to  be  about  1=  3  due  to  slipstream  effect. 
This  value  was  verified  experimentally  as  described  in  the  test  procedures 
section  and  found  to  be  a  reasonable  assumption.      The  lift  curve  slope  of 
the  tail  can  be  estimated  from  Reference  6  to  be,   a     =2.7/  rad. 


Cm,     =  -(2.7)(1.3)(.399)  =  -1.397 


A  lag  in  the  downwash  reaching  the  tail  may  be  found  from  the  emperically 
determined  formula:     (Reference  14) 

3/2 
CT      27  cos    '      A 
de 


Aa  =;7.3(.525  4  .475V736)AR-725(ll)-2e-385^ 


Where  y  is  the  vertical  distance  of  the  tail  quarter  chord  from  the  wing 
chord  line  (+  above). 
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d£ (4.  18)(27) 

d<J=  (57.31(4.  44C  725(2.  345)- 2e-°8°6 


da 


y  -   1,  0 
A=  0 

x  =  i.o 


Now  the  moment  change  with  angle  of  attack  change  may  be  computed  by: 


-       cCtW 


Cm      =  CT       x     +  C_        +Cm.  (1-1^-) 
ma  La     a  »f  xt  da 


C  =  (4.  18)(.  04)  +  .  0184  -  1.  397  (1   -  .  42) 


.6247 


«        (5.31)(.488)(.615) 


M       =  -12. 61 


There  is  also  a  moment  change  due  to  rate  of  change  of  angle  of  attack 
given  by: 
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lt  de 

M*    =  gky5    "T  d^  C^i 

vc   (5~)     c  t 


M      -  (32.2)(2,345)(.42)(-l,397) 
a  (84.45)(.488)(.615) 


M.    =  -1.746 
a 

The  principal  effect  of  pitch  is  manifested  in  a  damping  of  the 
moment  as  follows: 

u  M<*        -1-748 

Me  =  "^-=-~42- 

da 
M'Q  =  -4.  16 
The  elevator  control  power  may  be  expressed  by: 


T      C„ 


(St)      e    mit 

c" 

where  T  ,   the  control  effectiveness  parameter,    is  .6  from  Figure  5-33  of 
e 

Reference  7. 

(32.2)(.6)(-1.397) 
6  (.488)(5.31)(.615) 


16. 


Using  the  stability  derivatives  above,    the  following  determinant  of 
coefficients  is  obtained: 


AV 

La 

A8 

+  .236 

-6.2 

32.  2 

00903 

s  +  3.265 

-s 

00806 

1. 746s  +  12. 

61 

s(s+4.  16) 

From  this  the  characteristic  equation  for  the  rigid  aircraft  is: 

X4   +  9.  407X3   +  28.  4124Xa   +  7.  1315X  +  4.  514  =  0 
Solution  of  the  equation  yields  phugoid  and  short  period  characteristics  of: 


Phugoid 
C    =        .245 
uon    =        .414 
T|  =     7.  05  sec 


Short  Period 
C  =  -897 
tt>n  =  5.  128 

Tfc  =     .  151   sec 


Periods  15.7    sec  Period    =  1.22     sec 

ujc  =        .  401  cue  =  2.28 


Static  Stability 

The  neutral  point  and  maneuver  point  for  the  rigid  aircraft  can  be 
determined  from  the  stability  derivatives  once  they  have  been  converted 
to  non-dimensional  form. 
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I     -   X 

o       eg 

C 
m 
a 

a 

c 

m 

r    -  x 

m       eg 

C 
m 
a 

C 
mD9 

"CT. 

2 

C  =  -.6247 

m 
a 

C         =  .0342 


C  =  -.  3187 

mD6 


.29 


The  resulting  stick-fixed  neutral  and  maneuver  points  are: 

N     =  47.  5%  mac 
o 

N      =5  9.  9%  mac 


The  slope  of  the  neutral  point  curve  may  be  determined  by  differentiating 
d6 


j^—    by  x      . 
dC.  y     eg 
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d6  d6       dC  dC 

e    _         e  m    _         m 

dC~    "  dC~      dCT      "   dCT       C~ 


C 


'  6 


a  V 


m  m 

C      _(  Ct        "    2CT   * 
m.  L,  L 

6  a 


de 

CT    (x      -x      )  +  C™      +  C-.  (1-—)  C 

L,       eg       ac            "^z        mL.v        dc  m 

1       r        a        s                                ft  ^ 

C          L                               CT  "   2CT      J 

m_                                         L.  L, 


d6  / 

d(-^~-)/  _      1 


vdC 


L/dx  C  T     C^ 

eg  m  e      x 


With  the  elevator  effectiveness  estimated  as,   T     =  .  6,   the  slope  of  the 

e 

neutral  point  curve  is: 


d6 


-1. 191 


The  slope  of  the  maneuver  point  curve  may  likewise  be  deter - 
d&e 
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d6 

d6 

dC 

d6 

e 

e 

m 

dn 

dC 

dn 

dC 

dC 

m 

7H 

w 


W            m   _       W  m    _1_ 

C  dL      "  C  d~C.       qS 

m                           mB  J_i 
0                             6 


Cmtt{  +  Cmit  (1   "  5J) 

] 


The  slope  of  the  maneuver  point  curve  for  the  rigid  case  is 


d6 


d^dn-/dx        =-581 


Analog  Computer  and  Root  Locus 

Due  to  the  flexibility  of  this  aircraft,   the  vertical  bending  of  the 
rear  portion  of  the  fuselage  will  manifest  itself  in  a  change  in  the  hori- 
zontal stabilizer  incidence,    i  .     This  change  can  be  considered  as  a  change 
in  the  derivative  Cm.   .      This  effect  has  serious  aerodynamic  consequences 
because  all  the  pitching  moment  stability  derivatives,   with  the  exception 

of  speed  stability,    are  a  function  of  Cm,  . 

H 


20, 


The  effect  of  wing  bending  predominately  manifests  itself  in  a  change 
in  C^     .      In  this  investigation  a  varying  Ct       will  be  used. 

In  order  to  observe  the  result 
factored  out  of  the  moment  equation. 


V  M    Au  -  M  .  Ac.    -  M    Aa   +  A0     -  Mu  A0  -  Mc    A6     =0 
o     V  a  a  0  6        e 


t 

-  V  M_.Au  -  {  g,       a    —  ~  Cm.  J  A* 

°     V  VC    £*-)       c    dQ  H 


fer  fawxa  +  C^  +  C^  (1    -£)} 


+  Ae-{      S         _L  c      3a9-{ Wr  c      3  =  o 

vcT(^)     c         H  vcTF*)  H 


AG     =  -  V   M    Au  -  S         a    [a    x     +  C^      ]  Aa 

°     V  r*    "~  /_____.  w   a  a  f 


V-^Cm-   {(^!r-)A^  +  <1  -^)Aa+~A8.-T   A6     } 
ky         mit    v  V  da  do-  V  e      e 


This  leaves  the  moment  equation  in  a  suitable  form  for  both  root  locus  and 
analog  simulation. 
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The  root  locus  equation  for  varying  Cm.    was  derived  from  the 
K                 £ 
following  determinant  where  —  =  : a" 

2      CTr(^) 


D      -  g 
a       B 

V 
o 

L 


Q-  f  t 

4 


r    t  de  .,       de  .-,  , 


Upon  substituting  into  the  above  determinant  for  all  values  except  C         and 

a 
Cm.   ,   the  following  determinant  was  formed. 


S  +  .236  (.0736CL      -.3815)  +.3815 

a 

.761  s  +  .  781  CL  -  s 


.681  {.814  CT       =.3625  s(s  -  2.  977  Cm.  ) 

■Li  L 


-  C:        [1.251  s  +  11.69]} 
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This  produced  the  following  root  locus  equation  for  a  varying   Cm. 


4.238  Cm.  |[s3  +(2.  9943+.  5486CT      )s2  +(.  9406  +  .  0901  CT      )s+. 


8008] 


s4+(.236  +  .  781CL    )s3+(-.  0722=.6777CL    )s3+(=.  0855  -.  14CL  )s+(-.  1052  - 

.0292CL    ) 


Ct       was  not  substituted  in  because  the  root  locus  will  be  drawn  for 
a  range  of  CL     ls.     This  procedure  was  used  because  Cl      has  a  large  effect 
on  the  dynamics  of  the  aircraft.     Also  the  variable  C^      obtained  from 
Reference  3  could  not  be  depended  on  since  the  configurations  do  not  match 
exactly.     The  root  locus  will  be  drawn  for  a  range  of  possible  Cl    's, 
CT       =  3.  0,    4.  0,   and  5.  0  per  radian. 


=  3.0 


4.238  __.      [s3+  4.  6401s2  +  1.2110s  +  .8008] 

I     mitl 
__ _____ — __ » __    =   _i 

s4+  2.579s3-  2.  1053s2-  .  5055s  -  .  1925 

4.238|cm.  t(s  +  4.4065)(s  +  .  1168  ±j  .4099) 
(s  +  3.  1941  )(s  -  .  85  91  )(s  +  .  122  ±j  .2351) 
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4.23  8  C_   (s  +  5.  1887s  +  1.3  012s  +  .  8008) 
s4+  3.  36s3  -  2.  783s2  -  .  6455s  -  .2216 


4.23  8  Cm   (s  +  4.9588){s  +  .  1 149  ±  j  .  3852) 

L_v __, 

(s  +  4.  0166 )(s  -  .  8906 )(s  +.  117  ±j  .  2197) 


5.  0 


4.  238Cm      (s    +  5.  7373s**   +  1.  3913s  +  .  8008) 

I       Hi 

__ — _ — , — _ — _____  _   .j 

s4+  4.  141s3    -  3.4607s2  -  .  7855s  -  .2507 


4.  238  C_      (s  +  5.5112)(s  +  .  1131  ±j  .357) 

=  =1 

(s  +  4.  8268)(s  -  .  9143)(s  4  .  1143  ±  j  .  2092) 


The  root  locus  for  varying  C_.     for  the  three  Ct      !s  is  shown  in 
y     s      mit  '-'a 

Figure  16.     Figure  17  shows  the  stable  phugoid  portion  of  Figure  16  with 
an  expanded  scale.     Lines  of  constant  Cm.    are  shown  plotted  on  these 
curves.     Figure  17  also  shows  lines  of  constant  damping  ratio,  Q. 

The  analog  simulation  was  based  on  the  normal  lift  and  drag 
equations  plus  the  modified  moment  equation.      The  equations  used  in  the 
analog  simulation  are  shown  on  the  following  page. 
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\c  \r  i 

■  A9=   -J*{.V   M    Au  +-r-x   CT     Aa+^C  Ac*-  — — -  p]  dt 

J        o     V  2      a    L  2       m^,  2 


I  I 

p=  -L[K~  t^  A^  +K(1  -|^)Aa  +K-^  A0  +  KT  A6    ] 
V    da  da  V  e      e 


Au=-/{Dv-Tv)Au+^|fSICL    A,    --X.A.   +  ^-Ae}< 


Aa   =  -  E  LvAu  +       s  C       Aa    -  A0 

L    o  a- 


■  Ah  =  -  I'  {  V  A6  -  V    Aa  }  dt 
d         o  o 


Figure  13  shows  the  computer  set-up  used.     Table  II  is  a  tabulation  of  the 
pots  versus  the  quantity  they  set.     The  computer  set-up  was  made  in  such 
a  way  as  to  keep  Cl      and  Cm.     separate.     This  could  be  accomplished 
easily  for  Cmi    by  taking  advantage  of  the  previously  factored  equations. 
However,   the  Cl      could  not  be  so  easily  separated  since  it  enters  the 
lift,   drag,   and  moment  equations. 

The  computer  set-up  was  checked  by  comparing  the  analytic  solu- 
tion with  the  computer  solution.  This  Sanborn  tape  is  shown  in  Figure  14. 
As  can  be  seen  in  Figure  17,  the  computer  solution  plots  very  close  to  the 
analytic  solution. 
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Analog  computer  runs  were  made  for  Cl    's  of  3,  0,    4,  0  and  5.  0 

per  radian.     Cm,    was  varied  from  2.  0  down  to  the  Cmi    for  critical 

H  t 

damping.     Typical  Sanborn  tapes  are  shown  in  Figure  15. 

Agreement  between  the  root  locus  and  the  analog  computer  was 
checked  by  making  analog  runs  at  C^     !s  of  3.0,    4.  0,    and  5,  0  per  radian 
and  with  Cm-     =1.0.     The  cu    's  were  taken  directly  from  the  tapes  while 
the  C's  were  found  by  taking  amplitude  ratios  from  the  tapes  and  applying 
them  to  Figure  18.     These  values  are  shown  plotted  on  Figure  17. 

All  analog  runs  are  for  a  6     of  .  15  radians  introduced  as  a  pulse 
of  about  3  seconds  duration. 
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FLIGHT  TEST  PROCEDURES 


The  flight  test  phase  of  the  investigation  began  in  January  1964 
and  continued  through  early  May  1964  consisting  of  33  hours  of  flight  time, 
These  tests  were  conducted  at  the  air  facility  at  the  Forrestal  Research 
Center  of  Princeton  University.     This  phase  consisted  of  preflight  calibra- 
tions,   in-flight  calibrations,    stick  fixed  longitudinal  static  stability  flight 
test,    stick  fixed  longitudinal  dynamic   stability  flight  test,    flight  test  de- 
termination of  pitching  moment  due  to  tail  angle  of  attack  (Cm.   ),   and 

elevator  effectiveness  (T    ). 
e 

Preflight  Calibrations 

As  soon  as  the  test  instrumentation  was  installed  the  aircraft  was 
weighed  with  various  fuel  loads  to  determine  fuel  weight  and  moment  con- 
tributions.    The  aircraft  was  then  weighed  with  no  fuel  and  each  of  the 
three  pilots.     The  fuel  weight  and  moment  were  then  added  to  pilot  weight 
and  moment  to  arrive  at  center  of  gravity  locations  and  gross  weight 
versus  fuel  load  (Figure  19). 

It  was  found  that  three  pilots  weighing  170  lbs.  ,    190  lbs,   and 
210  lbs.  ,    respectively,    produced  center  of  gravity  excursions  from 
24.  10%  MAC  to  3  0.  00%  MAC.     An  attempt  was  made  to  extend  this  range 
of  center  of  gravity  locations  aft  by  adding  weights  to  the  tail;  however, 
this  introduced  nonlinearities  in  fuselage  bending  that  could  not  be  ac- 
counted for  in  the  static  stability  test.     It  was  decided  that  with  a  large 
selection  of  flight  test  points  the  static  stability  curves  could  be  ac- 
curately extrapolated  with  the  center  of  gravity  excursion  provided  by 
different  pilots. 

The  elevator  motion  was  calibrated  by  attaching  a  protractor  to 
the  elevator  and  recording  deflections  obtained  by  selecting  control 
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stick  positions.     This  calibration  was  found  to  be  different  for  7  psi    and 
5  psi    internal  pressures.     As  elevator  calibration  was  critical  for  static 
stability  tests,   the  elevator  calibration  was  done  every  time  flight  test 
data  was  obtained,     A  possible  explanation  for  the  change  of  elevator 
calibration  is  that  the  string  attached  to  the  potentiometer  could  change 
length  due  to  aging  and  ambient  weather  conditions.     The  use  of  wire  would 
eliminate  the  problem,   but  was  ruled  out  for  flight  safety  reasons.     A 
string  was  chosen  instead  of  a  wire  to  preclude  the  possible  jamming  of 
the  control  system  by  the  instrumentation. 

The  only  calibration  required  of  the  angle  of  attack  vane  was  to 
assure  that  its  nominal  position  fell  in  the  center  of  the  recorder  trace 
as  no  quantitative  results  were  required  of  it. 

The  recorder  accelerometer  was  simply  recorded  at  +1  and  -1  g 
in  order  to  check  the  flight  test  calibration  for  a  2  g    range. 

In  predicting  Cm,    for  the  analytical  development  it  was  necessary 
to  choose  a  value  for  tail  efficiency  (n  ).     It  was  felt  that  r\    for  this  air- 
craft  would  be  greater  than  1.  0  because  there  are  no  fuselage  effects  to 
dissipate  the  wake  other  than  that  obtained  by  the  downwash  effect  off  the 
wing  lowering  it  into  the  fuselage.      Considering  no  downwash  effect  the 
following  development  was  made: 


From  the 
momentum  theory: 


where  v  is  velocity  induced  by  disc  at  the  disc 
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(2)    V        =  V  +  2v  (downstream) 
ss 

Therefore  mass  flow  =  pA    (V  + 

and  from  (2): 

V        -  V 


Therefore: 


PAD 
mass  flow  =  — - —  (V        +  V) 


PAD 

T  =  -zr—  (V        +  V)(V        -  V) 


ss         PAD 

and  from  the  analytical  development  made  at  test  conditions,    i„  e.  , 
V  =  54  knots  IAS 


Drag  =  129.4  lbs  =  Thrust 


Therefore: 


v       =   /      (2X129.4)  f 

ss     N/(0.0024)(12.02)         l°    'D) 


V        =  127  fps 
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Therefore  from  (2): 

2v  =  42.5  fps 

Therefore  if  no  downwash  effect  is  considered  and  accepting  the 
assumptions  that  went  into  the  momentum  theory  (i.  e.  ,    constant  v  across 
disc,    etc.  ),   the  difference  between  the  velocity  in  the  wake  and  the  free 
stream  should  be  42.5  fps.     Assuming  this  the  case,   the  engine  was  run 
statically  and  power  adjusted  to  produce  maximum  wake  velocities  of 
40,    45,    and  50  fps.     These  surveys  were  made  with  a  velometer,   which 
is  a  mass  flow  device.     At  these  power  settings  wake  surveys  were  made 
at  the  tail  with  the  velometer  and  average  tail  velocities  plotted  versus 
wake  velocity  on  Figure  20. 

Entering  Figure  20  with  wake  velocity  of  42.  5  fps  an  average  tail 
velocity  of  31.2  fps  was  obtained.     This  velocity  is  then  the  prop  wake 
contribution  to  tail  velocity  at  the  flight  condition.     From  this  tail  effi- 
ciency was  calculated  as  follows: 


Ht        (84.5  +  31.  2f 


<1  (84.  5)2 

Tit  =  1.88 

This  figure  is  most  idealized  in  that  many  assumptions  were 
necessary  to  obtain  it;  however,    it  does  justify  values  of  r\    greater  than 
1.0.     If  the  downwash  effect  is  considered  an  T)     =  1.  3  is  not  considered 
to  be  out  of  order.     Consequently,    this  value  was  used  in  the  analytical 
development. 

It  was  also  necessary  to  estimate  the  change  in  Cm.     due  to  fuse- 
lage bending.     The  aircraft  was  suspended  and  a  25  lb  weight  placed  on 


the  tail.     At  internal  pressures  of  5  psig  and  7  psig  the  deflection  of  the 
tail  relative  to  the  wing  root  chord  was  found  to  be  1.8  degrees  for 
5  psig  and  0.  9  degrees  for  7  psig  internal  pressure.      Therefore: 


[) 


La 

>.rf/'7  0    IF  =TT8des/lb-f 


r) 


La 

~)        =~<ieg/lb.f 

D  .  U 


Then  for  any  change  in  force  on  tail  AF: 


AF 
iart'7.0=    27^8  d6g 


AV5.0  x  TT9deg 


due  to  bending 

Knowing  to     and  q(8.  46  psf  at  S,  L.)a  change  in  a     of  1  degree  will  result 
in  a  change  in  force  on  the  rigid  tail  of  8.  4  lb.  f . 
Therefore: 


AF  =  8.  4  La 
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Now,    the  change  in  angle  of  attack  of  the  tail  due  to  the  flexible  fuselage 
may  be  expressed  as: 


La     =  La  -  La 

t  t    .    .  .  t.  , 

rigid  due  to  bending 


AF  A 

La      r  —-   -  La 

t        K  t    .    .. 

rigid 


8.4   . 
La     t  — —  La     =  La 

t  K  t  t    .    .. 

rigid 


Solving  for  La   : 


rigid 


Since  C  varies  directly  with  La     and  substituting  for  values  of  K 

found  previously: 

C         )  .  0.  768  Cm     ) 

1    1.0  rigid 


t)      =0-623cmi) 

4   5.0  frigid 
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From  analytical  development  Cm  =  1.  397 

Xt    r.gid 


Therefore: 


C         )         =  1.075 

x    7.  0 


Cmi)  -0.871 

^5.0 


These  values  will  be  compared  with  other  methods  of  C  deter  = 

mination  in  the  results  section  of  the  report. 

Flight  Calibrations 

The  in-flight  calibrations  consisted  of  airspeed  and  normal  ac- 
celeration calibrations.     The  cockpit  airspeed  indicator  calibration  was 
obtained  by  flying  the  aircraft  at  approximately  20  feet  above  the  runway 
at  indicated  airspeeds  of  45,    50,    55,    603   and  65  knots.     These  runs  were 
made  on  reciprocal  runway  headings  in  order  to  compensate  for  wind. 
Knowing  the  runway  length  the  average  time  for  each  run  enabled  a  true 
airspeed  calculation.     This  true  airspeed  was  then  corrected  for  tempera- 
ture and  pressure  altitude  to  obtain  calibrated  or  corrected  indicated  air- 
speed.    The  results  of  the  calibration  are  shown  in  Figure  21.     Calibra- 
tions were  run  on  January  3rd,    January  7th,   and  April  25th  with  agree- 
ment on  all  three  dates.      The  recorder  airspeed  indicator  was  calibrated 
against  the  cockpit  indicator.     However,   as  mentioned  previously,    the 
static  source  for  the  recorder  airspeed  indicator  was  a  venturi  mounted 
on  the  angle  of  attack  boom  mount.     The  aircraft  pilot  static  tube  was 
used  as  the  dynamic  pressure  source  for  each  instrument. 

The  recorder  accelerometer  calibration  was  obtained  by  readings 
made  of  two  cockpit  accelerometers  in  steady  turns.     The  aircraft  was 
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equipped  with  a  standard  cockpit  accelerometer  and  a  glass  tube  accele - 
rometer  was  also  added  as  a  check.     The  two  instruments  were  in  agree- 
ment and  the  calibration  of  the  recorder  accelerometer  did  not  change 
throughout  the  flight  test. 

Flight  Test  of  Stick  Fixed  Static  Stability 

The  static  stability  tests  were  conducted  at  a  power  setting  to  give 
54  knot  IAS  in  level  flight.      (This  was  the  condition  used  in  the  dynamic 
analysis,  )     This  power  setting  was  held  constant  for  all  flight  tests, 

An  accurate  record  of  fuel  on-ha-rtd  was  kept  in  order  to  know 
weight  and  center  of  gravity  location  for  calculation  of  lift  coefficients 
and  correcting  normal  acceleration  to  a  standard  gross  weight.  Fuel 
was  measured  at  the  beginning  and  end  of  each  flight.  A  consumption 
rate  of  4  gallons  (24  lbs)  per  hour  was  obtained  and  this  rate  was  used 
to  calculate  fuel  on  board  for  each  recorded  event. 

Determination  of  Stick  Fixed  Neutral  Points 


Neutral  points  were  obtained  for  internal  pressures  of  7.  0  and 
5.  0  psig.     The  flight  procedure  was  to  set  the  power  at  54  knots  IAS  in 
level  flight.     Elevator  deflection  and  airspeed  were  then  recorded  at 
indicated  airspeeds  of  45,    50,    55,   and  60  knots  after  a  steady  state 
flight  condition  had  been  reached.      These  runs  were  repeated  by  each 
pilot  to  obtain  data  at  three  different  center  of  gravity  locations. 

For  the  7  psig  runs  an  additional  point  was  taken  at  65  knots  IAS, 
however,   at  5  psig  an  undamped  high  frequency  oscillation  was  en- 
countered at  63  knots  IAS.     An  explanation  of  this  motion  is  presented 
in  the  Discussion  of  Results.     Figure  22  shows  a  typical  trace  for  ele  = 
vator  versus  airspeed  and  also  a  trace  of  the  oscillation. 
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Elevator  deflection  was  then  plotted  versus  lift  coefficient  for 
each  center  of  gravity  location.      These  plots  were  linear  and   the  slopes, 


point.     Figures  23  through  29  show  the  plots  of  6     versus  C      and 

e  L 

Figures  30,  -3-r7  and  32  show  plots  of  6    /C      versus  X       for  5  psig,    7  psig 
and  rigid  calculations. 


Determination  of  Stick  Fixed  Maneuver  Points 

Maneuver  points  were  also  obtained  for  internal  pressures  of 
5.  0  psig  and  7.  0  psig.     For  these  test  power  was  set  at  54  knots  IAS  in 
level  flight.      The  aircraft  was  then  slowed  to  45  knots  IAS  and  pushed 
over  to  obtain  5  7  knots  IAS.     A  pull  up  was  then  initiated  with  a  constant 
elevator  position.     During  the  pull  up  airspeed,    normal  accelerations,   and 
elevator  position  were  recorded.      The  tests  were  repeated  for  each  cen- 
ter of  gravity  location.     At  the  point  on  the  trace  where  the  airspeed  was 
54  knots  IAS,   acceleration  and  elevator  positions  were  read  and  plotted 
on  Figures  23  through  29.     Figure  22  shows  a  typical  trace  of  6    ,  V„   and 

n.     The  normal  acceleration  was  corrected  to  a  standard  gross  weight  of 

S7S~ 

§£tT  lbs  in  order  to  match  the  flight  condition  used  for  the  dynamic  calcu- 

lations  (i.  e.  ,    52-0  lbs  gross  weight,    center  of  gravity  at  29.  00%  MAC, 

and  54.  0  knots  IAS).      The  plots  of  6     versus  n  were  linear  up  to  approxi- 

e 

mately2g's.  At. load  factors  above  this  the  fuselage  bending  became  non 

linear.     However,    since  the  dynamic  motion  was  later  found  to  be  in  the 

linear  region,    these  linear  slopes  were  used  for  plotting  d6    /  dn  versus 

X       to  obtain  maneuver  points.     Figures. 3-©t   31,   and  32  show  plots  of 

d6    /  dn  versus  X        for  5  psig,    7  psig,   and  rigid  calculations, 
e  eg 
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Determination  of  Longitudinal  Stick  Fixed  Dynamic  Stability- 
Dynamic  tests  were  conducted  at  5  psig  and  7  psig  internal  pres- 
sures.    For  the  dynamic  test  the  aircraft  was  loaded  so  as  to  obtain  a  gross 

sr/j 

weight  of  &2rQ  lbs  and  a  center  of  gravity  location  of  29.  00%  MAC  in  order  to 
match  the  condition  of  the  analytical  development. 

The  aircraft  was  excited  with  a  pulse  elevator  input,     It  was  found 
that  some  method  of  returning  the  pulse  to  the  original  amplitude  and  hold- 
ing it  there  was  necessary.     In  order  to  achieve  this  a  chord  was  attached 
to  the  pilot's  lap  belt  and  to  a  ring  that  could  be  slipped  over  the  control 
stick.     The  length  of  the  chord  was  adjusted  to  maintain  54  knots  IAS  in 
level  flight.     After  each  pull  up  the  control  stick  was  held  against  the 
chord  to  obtain  a  constant  elevator  deflection.     This  sort  of  "Rube  Goldberg" 
arrangement  is  typical  of  what  must  be  resorted  to  when  an  airframe  will 
not  accept  screws  or  nuts  and  bolts  graciously. 

During  the  dynamic  test  elevator  deflection,   airspeed,   angle  of 
attack,    pitch  rate,   and  acceleration  were  recorded,     Figure  33  is  a  typical 
recorder  trace  showing  the  dynamic  response. 

In  addition9    during  some  of  the  tests  the  gun  camera  recorded 
fuselage  deflection  during  the  phugoid  motion,     A  flutter  of  the  elevator 
was  used  to  correlate  events  on  the  camera  and  the  recorder  trace.    From 
the  recorder  trace  phugoid  damping  ratio  (C)  was  computed  by  Figure  18. 
Damping  ratio  and  frequency  of  phugoid  oscillation  (cu    )  were  then  plotted 
on  Figure  17  which  is  a  root  locus  plot  obtained  from  the  analytical 
development. 

The  short  period  oscillation  was  not  investigated.     The  analytical 
development  indicated  at  short  period  time  constant  of  the  order  of  0,  2 
seconds.     The  angle  of  attack  vane  would  not  respond  at  this  rate  and  from 
matching  elevator  input  and  angle  of  attack  response  it  is  impossible  to 
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differentiate  between  aircraft  short  period  and  angle  of  attack  vane  response, 
In  any  event,    the  highly  damped  short  period  mode  was  not  considered  in- 
teresting enough  to  warrent  a  more  intense  investigation. 

Movie  shots  of  fuselage  bending  were  projected  on  a  grid  screen. 
From  these  observations  it  was  possible  to  determine  fuselage  bending,    in 
inches,    at  the  tail  and  where  in  the  motion  the  bending  occurred,     From  the 
preflight  test  calibrations  it  was  determined  that  the  maximum  bending  for 
the  5  psig  runs  was  of  the  order  of  1 ,  5  degrees,   and  1.  0  degrees  at  7,  0 
psig.      This  order  of  bending  agrees  with  the  magnitudes  found  in  the  pre- 
flight test  determination  of  Cm.    by  statically  bending  the  fuselage.     Fig- 

xt 
ure  34  shows  three  typical  frames  from  the  movie  shots. 

Determination  of  C™,    and  T 


In  order  to  further  verify  the  analytical  development  it  was  neces- 
sary to  find  the  aircraft  pitching  moment  due  to  tail  incidence  change  (Cm.   ) 

it 
and  elevator  effectiveness  (T  ). 
e 

To  calculate  Cm.   ,   flight  test  data  from  Figures  23  through  29 

was  used  to  obtain  values  of  6     for  center  of  gravity  locations  at  a  constant 

C      of  0.  5  (this  was  the  value  used  for  the  rigid  analysis).     In  this  manner 

a  slope  AX      /A6     was  obtained  from  Figure  35, 
eg  e 

From  previous  flight  tests  it  was  determined  that  the  aircraft  tail 

incidence  was  too  high  in  that  up  elevator  was  necessary  for  almost  all 

flight  conditions  (note  Figure  2),     Therefore,   after  all  other  flight  tests 

were  completed  the  tail  incidence  was  lowered  5  degrees.     The  aircraft 

was  then  flown  at  C      of  0,  5  at  three  center  of  gravity  locations  and  elevator 

position  recorded.     As  AX      /A6      should  not  change  with  tail  incidence  a 
eg         e 

line  was  faired   through  these  points  with  the  AX      /  A5     slope  found  pre- 
viously. 


From  Reference  13: 
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m6  L 


A6 


(IfC^         =  cL       ) 

wing  4a4i  £/<, 


Therefore: 


and: 


aL       r        A6 


r       -  t     -  5   de§° 

e       A6  A6 


From  Figure  35  T     was  found  to  be  0.  55.     Using  this  value  of  T    ,    the  slopes 

AX      /  A6     from  Figure  35,   and  equation  1: 
eg         e  s 

C  =  =1.0 

*5.0 


C  -  =1.23 

*7.0 


These  values  will  be  discussed  in  the  results  section  of  the  report. 
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Other  Aspects  of  Flight  Test 

Some  of  the  other  unique  features  of  the  aircraft  effected  flight 
test  procedures.      The  semi-open  cockpit  induced  a  chill  factor  that  had  a 
decided  effect  on  pilot  response  during  the  late  winter  months.     Also,   the 
open  cockpit  proved  a  problem  for  voice  transmission.     This  was  solved 
by  wearing  an  oxygen  mask  with  valves  removed,      This  procedure  cut  the 
background  noise  to  a  very  low  level. 

Due  to  the  aft  mounted  propeller  any  loose  gear  in  the  cockpit 
could  easily  reach  the  propeller  arc  and  be  responsible  for  a  fabric  punc- 
ture.    Thus,    everything  in  the  cockpit  was  of  necessity  tied  down. 

The  low  aerodynamic  mass  of  the  aircraft  made  it  very  suscep- 
tible to  light  gust  ,   and  to  obtain  usable  data  extremely  smooth  air  was 

33 
required.     Figure«2-5r*  shows  a  noise  trace  on  a  day  that  would  normally 

be  considered  quite  smooth. 

Ground  handling  of  the  aircraft  was  reminiscent  of  sea  plane 
operation  and  was  a  three  man  job,  The  ability  of  the  aircraft  to  use 
sod  runway  greatly  simplified  the  flight  test,    however. 

In  all  respects  the  three  pilots  concerned  were  quite  pleased  to 
be  able  to  obtain  a  small  taste  of  open  cockpit  and  the  scream  of  wires 
that  was  so  typical  of  early  aviation, 
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DISCUSSION  OF  RESULTS 

Definition  of  Problem 

The  original  purpose  of  this  investigation  was  to  determine  the 
applicability  of  standard  analytical  techniques  on  predicting  static  and 
dynamic  stability  of  an  inflatable  airplane.     It  was  assumed  that  the 
major  factor  influencing  any  deviation  from  rigid  analysis  for  this  air- 
craft would  be  fuselage  bending.     A  secondary  effect  would  be  wing  tor- 
sion and  bending  affecting  lift  curve  slope  (C,      )  with  changes  with 
dynamic  pressure  (q).     The  change  in  Cj_,      with  q  was  obtained  from 
Reference  4  which  is  a  report  on  full  scale  wind  tunnel  test  of  the  air- 
craft.    As  this  data  was  available  and  accounts  for  flexibility  effects  on 
the  wing  and  tail,    it  was  decided  to  concentrate  on  the  effects  of  fuselage 
bending. 

The  result  of  research  was  that  there  are  many  different  opinions 

as  to  the  effect  of  fuselage  bending  on  such  things  as  movement  of  stick 

fixed  neutral  and  maneuver  points  (N     and  N     ).     It  was  then  determined 

o  m 

that  a  further  investigation  of  the  effects  of  fuselage  bending  on  longi- 
tudinal stability  might  be  profitable.     This  aircraft  would  seem  to  be  an 
ideal  vehicle  for  this  type  of  investigation  in  that  fuselage  stiffness  could 
be  varied  by  changing  internal  pressure. 

All  references  agreed  that  fuselage  bending  was  an  overall  de- 
stabilizing effect,   and  were  in  general  agreement  that  N      would  move 

forward,    except  for  Reference  9.     Reference  9  stated  that  N      would 

r  m 

move  aft  due  to  a  change  in  aircraft  center  of  lift.     The  movement  of 

N     was  not  in  agreement,    however,    in  that  there  was  difference  of 
o  s 

opinion  in  how  far  N     should  move  forward  relative  to  N 

o  m 

It  was  anticipated  that  the  change  in  N      and  N     between  internal 
pressures  of  5  psig  and  7  psig  would  answer  some  of  these  questions. 
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It  was  considered  that  a  valid  measure  of  fuselage  bending  would 
ages  in  Cm.   , 
fuselage  flexibility. 


be  changes  in  Cm.   ,   and  it  was  expected  that  Cm.    would  be  reduced  by 


Predicted  Results 

The  rigid  analysis  was  made  using  classical  techniques  of 
References  2  and  7.     An  effort  was  made  to  determine  tail  efficiency  (r|  ) 
by  a  wake  survey  mentioned  in  the  preflight  test.      This  yielded  a  value 
of  r]    of  1.  3  that  resulted  in  a  Cm.    for  the  rigid  aircraft  of  -1.  397.     A 
full  scale  wind  tunnel  value  of  Cj^     was  chosen  at  the  test  q  of  8.  46  psf. 

The  results  of  this  analysis  were  a  predicted  N     of  47.  5  percent 

MAC,    N       of  59.  9  percent  MAC,    phugoid  period  of  15.  7  seconds,   and  a 
m 

damping  ratio  of  0.245. 

For  the  semi-rigid  analysis  Cm.     and  Cj^      were  factored  out  of 
the  equations  of  motion  and  a  root  locus  plotted  showing  effects  of 
changing  Cm.    and  Cj^      (Figures  16  and  17).     The  excursion  of  C^,     was 
expected  to  be  from  3.  0  to  5.  0  and  values  of  3.0,    4.  0,   and  5.  0  per  radian 
were  used  for  the  root  locus  plot. 

A  static  bending  test  of  the  fuselage  was  made  at  internal  pres- 
sures of  5.  0  psig  and  7.  0  psig.      (This  procedure  may  be  found  in  the 
flight  test  procedures  section  of  the  report.  )    The  results  of  this  test 


Flight  Test  Results 

Although  the  flight  test  determination  of  Cm.     and  T     were  done 
at  the  end  of  the  program  a  flight  test  value  of  T    '^necessary  to  corre- 
late other  flight  test  results  with  predictions.     From  Figure  35,    T     was 
found  to  be  0.  55.    (The  method  used  may  be  found  in  the  flight  test  sec- 
tion of  the  report.  )     This  value  of  T     compared  favorably  with  the  pre- 

e 

dieted  value  of  0.  6  and  is  considered  valid. 
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data  this  result  can  only  be  considered  as  a  qualitative  indication  in  the 
movement  of  N     with  fuselage  bending,     A  possible  explanation  for  the 
scatter  in  the  neutral  point  data  may  be  due  to  a  change  in  position  of 
the  elevator  potentiometer  string  tie  down  point.     As  velocity  is 
changed,    the  leading  edge  of  the  horizontal  stabilizer  may  deform  and 
twist,    consequently  displacing  the  tie  down  point.     Since  the  neutral 
point  tests  were  made  over  a  speed  range  of  45  to  65  knots,    the  tie 
down  point  would  not  remain  in  the  same  relative  position,    thereby 
making  the  elevator  readings  differ  from  the  exact  value  in  an  unpre- 
dictable manner.      It  can  be  concluded  that  N     moves  forward  with 

o 

fuselage  bending;  however,    the  use  of  this  data  to  conclude  the  rela- 
tive movement  of  NT     and  N      with  flexibility  is  questionable, 
o  m 

Although  the  slope  of  the  6    /C      versus  X        curve  is  not  con- 
e        L  eg 

sidered  valid  in  finding  Cm.   ,    for  a  non-rigid  aircraft,    values  of  Cm, 

were  calculated  from  Figure  30.      Cm.  was  found  to  be   -1.  33  and 

Cm.  -1.  54,    using  the  flight  test  value  of  T    .      Using  the  data  of 

the  6    /  CT    curves,    Figures  23  through  29,    as  plotted  on  Figure  35, 
e        L 

C^,.  was  found  to  be   -1.0  and  C-^.  -1.23.      This  procedure 

H  5.  0  H   7,  0 

is  described  in  the  flight  test  section  of  the  report  and  is  considered 

a  much  more  valid  method  of  determining  Cm.     for  a  non-rigid  air- 
craft than  using  the  slope  of  the  6    /  C      versus  X        curve.      These 
s  e        L  eg 

values  are  in  excellent  agreement  with  those  predicted  by  the  static 

bending  test. 

The  flight  test  results  for  N      were  found  to  have  much  less 
&  m 

scatter  and  are  considered  to  be  of  significance.      Nm  5^  q  was  found 

to  be  46.  0%  MAC  and  N      _    n    49.  2%  MAC.      These  results   show  a 
m  7 .  i» 

definite  destabilizing  trend  with  bending.      The  slopes  of  the  6    /n 
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versus  X        curves  of  Figure  31  yielded  values  of  C^  of  -1.  07  and 

Cg  mit5.0 

Cmi  of  -1-40.     These  values  are  in  good  agreement  with  magnitudes 

predicted  by  the  static  bending  test. 

Phugoid  periods  and  damping  ratios  obtained  from  rigid  analysis 
and  flight  test  data  were  plotted  on  the  root  locus  diagram,    Figure  17. 
The  results  of  these  tests  showed  a  definite  trend  towards  instability 
with  increasing  fuselage  flexibility.      Cmi  was  found  to  be   -0„  8  and 

Cm.  -1.  25  by  the  mean  position  of  the  plots  of  flight  test  periods 

and  damping  ratios  at  the  two  inflation  pressures. 

The  short  period  mode  of  dynamic  motion  could  not  be  investi- 
gated with  the  instrumentation  installed.      It's  high  damping  ratio  of  0.  89 
and  low  time  constant  of  0.2  seconds  precluded  its  usefulness  for  these 
tests.      The  reason  for  the  high  damping  ratio  is  the  very  low  aerodynamic 
mass  of  this  aircraft. 

Table  III  is  a  comparison  of  Cm.     values  found  from  static  bending 
tests  and  four  flight  test  methods.      This  comparison  shows  close  agree- 
ment by  all  methods  and  a  definite  destabilizing  trend  with  fuselage  bend- 


used  in  the  rigid  analysis.      It  can  then  be  stated  with  some  certainty 
that  the  rigid  analysis  was  valid.      This  being  the  case,    N     moved  for- 
ward from  the  rigid  value  to  the  5  psig  value  7,  5%  MAC  while  N       moved 

»  o  m 

forward  13.  9%  MAC.     While  the  movement  of  N     between  5  psig  and 

7  psig  is  not  well  defined,    its  movement  from  a  rigid  value  to  a  flexible 

value  is  considered  to  be  of  such  validity  as  to  warrant  mention. 

It  may  therefore  be  stated  that  for  this  aircraft  N       and  N     are 

m  o 

both  destabilized  by  fuselage  flexibility  and  it  is  further  indicated  that 

N       moves  further  forward  than  N    . 
m  o 

The  high  frequency  aircraft  oscillation  mentioned  in  the  flight 
test  procedures  was  determined  to  be  caused  by  elevator  flutter.      This 
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flutter  was  found  to  become  excited  at  63  knots  IAS  with  5  psig  internal 
airframe  pressure. 

A  movie  supplement  to  Reference  4  relates  a  similar  motion 
occurring  at  7  psig  pressure  at  a  dynamic  pressure  of  20  psf  in  the 
full  scale  wind  tunnel.     The  motion  was  felt  to  be  excited  by  elevator 
flutter  in  that  fixing  the  elevator  to  the  horizontal  stabilizer  stopped 
the  motion. 

Figure  22  shows  a  recording  of  the  motion  during  flight  test  at 
5  psig  internal  pressure.     From  the  recording  it  was  found  that  the  ele- 
vator oscillation  was  from  amplitudes  of  +  2  degrees  to  +  13  degrees  at 
a  frequency  of  5  cps.     It  is  felt  that  this  is  the  same  motion  described 
in  the  wind  tunnel  test. 

In  the  wind  tunnel  analysis  the  cause  of  the  flutter  was  attributed 
to  control  looseness.      This  assumption  is  felt  to  be  valid  in  that  the 
elevator  can  be  moved  approximately  5  degrees  with  the  stick  fixed  at 
7  psig  internal  pressure.     At  a  lower  pressure  the  elevator  is,    of 
course,    more  loose.     It  is  felt  that,   this  being  the  case,    the  elevator 
flutter  could  occur  at  a  dynamic  pressure  of  12  psf  (63  knots  IAS),   and 
that  the  motion  observed  in  flight  test  is  the  same  motion  previously 
seen  in  the  wind  tunnel. 
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CONCLUSIONS 

It  is  concluded  that  rigid  analysis  techniques  will  predict  longi- 
tudinal stick  fixed  dynamic  stability  of  the  YAO-3G  inflatable  airplane 
to  a  good  degree  of  accuracy.     Static  stability  predictions  are  not  as 
accurate  but  do  indicate  typical  values.     The  aircraft  when  operated  at 
recommended  internal  pressures  is  of  sufficient  rigidity  as  not  to 
deviate  far  from  rigid  analysis.     However,    if  rigid  values  of  Cm.    are 
corrected  by  factors  obtained  from  static  bending  tests,   accurate  pre- 
dictions may  be  obtained  for  both  dynamic  and  static  stick  fixed  sta- 
bility. 

The  effect  of  fuselage  flexibility  of  this  aircraft  was  destabiliz- 
ing to  both  dynamic  and  static  stability.     The  phugoid  mode  was  ren- 
dered less  stable  by  increased  flexibility.     Both  stick  fixed  neutral 
and  maneuver  points  moved  forward;  however,    it  was  indicated  that 
N       is  destabilized  to  a  larger  extent  than  N   . 
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TABLE  I 
Aircraft  Specifications 
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General • 


Span 
Length 
Height 
Wing  area 

Wing  loading 
Thrust  loading 
Moment  of  Inertia  (Iy) 


Weights: 


Engine: 


Weight  empty 

Useful  load  =  1  pilot 

Fuel 

Gross  take-off 

Available  sea  level  static  thrust 
at  combat  rating 

Propeller  diameter 

Height  of  thrust  line  above  CG 

Minimum  fuel  consumption  at 
sea  level  cruise 


23 

,2     ft 

22 

,  0     ft 

7, 

,2     ft 

124, 

,3     ft2 

4, 

,42  #/ft2 

1, 

,69  #/ft2 

277. 

2     slug  ft2 

315 

lbs 

194 

lbs 

10 

gal      61  lb 

570 

lbs 

220  lbs 

at  3750  RPM 

47 

in 

3. 

14  ft 

20. 

2  #/hr 
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TABLE  I  (Continued) 

Wing: 

Airfoil  section 

Area   (S) 

Span  (b) 

Aspect  ratio  (AR) 

Taper  ratio  (X) 

Chord  (c) 

Ailerons: 

Chord  c    /  c° 
AreaSa-f-S- 
Horizontal  Tail: 

Area  (S  ) 

Span  (bt) 

Chord     (c    )  , 

st  \ 

Tail  volume  (V    = — ) 

t      S     c 

Elevators: 

Chord  (c    )  1.33     ft 

e 

Area  8.86     ft2 

Elevator  travel  3  0  up  to  25  down 

Vertical  Tail: 

Area    (S    )  15.31     ft2 

v  -v 

Span  (b    )  3.  92     ft 

v 

Chord  (c    )  4.416  ft 

v 

Rudder: 

Mean  chord  1.  66     ft 

Area  5.0      ft2 


0015 

124. 

3     ft2 

23. 

2     ft 

4. 

44 

1 

5. 

31     ft 

20. 

4     % 

6. 

41  ft2 

21. 

05  ft2 

6. 

66  ft 

3. 

167  ft 

399 

TABLE  II 
Analog  Computer  Potentiometer  Values 
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Pot 

Quantity 

Value 

Pot  Setting 

1 

6 
e 

K    ^ 

.  1 

.  1 

2 

5.  954 

.  5954 

3 

KTe 

24.  18 

.2418 

4 

«<i-£> 

23.  38 

.2338 

5 

o 

2,  502 

.2502 

6 

Cmit/  2 

CmH/2 

CmH/2 

7 

V   M„ 
o     V 

.681 

.681 

8 

K  Cm 

,3625 

.  3625 

9 

K    - 

T"  x   a 
2       aw 

. 806  CT 

.  0806  CL 

10 

D     -  T 
V         V 

.  236 

,236 

11 

V 
o 

.  3815 

.  3815 

12 

V 
o 

.3815 

.  3815 

13 

e  IT  AR  aw 

. 0736  CT 

. 0736  CT 

14 

V 
o 

84.45 

.8445 

15 

c7V"aw 

L      o 

.  781  CT 

. 0781  CT 

16 

Lv 

,761 

.761 

17 

V 
o 

84.45 

.  8445 

18 

V 
o 

84.45 

.8445 

Voltage -quantity  ratios 

Aa,A9,  A6    :   1  v.    =  .  01   rad  Au,    Av:   1  v.    =  1  ft/sec 

'       *       e 
LatbS  :  1  v.    =  .  01   rad/  sec  Ah:   1  v.    =  1  ft 


TABLE  III 


Cmit  Values 


50. 


Rigid 

7. 0  psi 

5. 0  psi 

Rigid  Analysis  using 
Classical  Techniques 

-1. 397 

- 

- 

Static  Bending  Test 

(with  factors  applied 
to  rigid  values) 

- 

-1.  08 

-    .  87 

Static  Flight  Tests 
a)  gradient  of 

d6e             CL-1 
dn        TeCm.t 

Te         =  .55      CT          =.5 
eexp                     -hi=l 

Figure  31 

- 

-1.4 

-1.  07 

b)  gradient  of 

d6e                 1 

- 

-1.54 

-1.33 

dCL       TeC-it 

Figure  30 

C         A6 

c)C          =_L_         e 

it        T        Ax 

e           eg 

Figure  35 

- 

-1.23 

-1.  o 

Dynamic  Flight  Tests 
Figure  17 

- 

-1.25 

-    .  8 

Average 

- 

-1.30 

-1.  01 
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Figure  33 .  SFIM  Recorder  Tapes 


Level  Flight  Condition 


Deflection  Due  To 
Up  Load 


Deflection  Due  To 
Down  Load 


Fig.    34.    Selected  Prints   From   Movie  Camera  Showing  Fuselage  Bending. 
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